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EFFECTS OF PREHEATING ON X-RAY SCALING RELATIONS IN GALAXY CLUSTERS 
John J. Bialek 1 , August E. Evrard 1,2 & Joseph J. Mohr 3,4,5 
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ABSTRACT 

The failure of purely gravitational and gas dynamical models of X-ray cluster formation to reproduce 
basic observed properties of the local cluster population suggests the need for one or more additional 
physical processes operating on the intracluster medium (ICM). We present results from 84 moderate 
resolution gas dynamic simulations designed to systematically investigate the effects of preheating — an 
early elevated ICM adiabat on the resultant, local X-ray size-temperature, luminosity-temperature 
and ICM mass-temperature relations. Seven sets of twelve simulations are performed for a ACDM 
cosmology, each set characterized by a different initial entropy level Si. The slopes of the observable 
relations steepen monotonically as Si is increased. Observed slopes for all three relations are reproduced 
by models with Si <G 55 — 150 keV cm 2 , levels that compare favorably to empirical determinations of 
core ICM entropy by Lloyd-Davies, Ponman & Cannon. The redshift evolution for the case of a locally 
successful model with Si — 106 keV cm 2 is presented. At temperatures kT > 3 keV, little or no evolution 
in physical isophotal sizes or bolometric luminosities is expected to z < 1. The ICM mass and total 
masses at fixed T are lower at higher z as expected from the virial theorm. ICM mass fractions show a 
mild T dependence. Clusters with T < 3 keV contain ICM mass fractions depressed by modest amounts 
(< 25%) below the cosmic mean baryon fraction f2(,/fi m ; hot clusters subject to preheating remain good 
tracers of the cosmic mix of clustered mass components at redshifts z < 1. 

Subject headings: clusters: general — intergalactic medium — cosmology 



1. INTRODUCTION 

The hot, X-ray emitting plasma known as the intraclus- 
ter medium (ICM) represents a part of the baryonic mat- 
ter of the universe that is not associated with individual 
galaxies but remains trapped in the deeper gravitational 
potential of galaxy clusters. Models of cluster formation 
in which the intergalactic gas simply falls into the dark 
matter-dominated gravitational well (so-called infall mod- 
els) fail to reproduce all the structural properties of the lo- 
cal cluster population (e.g. Evrard & Henry 1991; Navarro, 
Frenk & White 1995; Mohr & Evrard 1997; Bryan & Nor- 
man 1998). There appears to be additional physics driving 
ICM evolution. 

As early ROSAT and Einstein data emerged, several 
researchers proposed that the missing element is the exis- 
tence of a high entropy intergalactic gas prior to a cluster's 
collapse (David et al. 1991; Evrard & Henry 1991; Kaiser 
1991; White 1991). The entropy floor preduced by the 
preheated gas breaks the self-similarity between the dark 
matter and ICM on different mass scales. This occurs be- 
cause the equivalent thermal energy of the entropy floor 
corresponds to a larger fraction of the virial temperature 
in lower mass clusters. 

A physical scenario consistent with this picture is one 
in which the bulk of star formation in a proto-cluster re- 
gion occurs early in its formation history. Heating from 
supernova-driven galactic winds and AGN activity asso- 
ciated with the rapid star formation ultimately comes to 



exceed the cooling rate of ambient material, resulting in a 
period of net local heating. In the low-density limit, for 
which subsequent cooling is unimportant, this feedback de- 
fines an initial adiabat Si that is conserved in the absence 
of shocks. Later shock heating can raise the adiabat above 
this "entropy floor," but the entropy cannot decrease be- 
low this level unless cooling becomes important. 

Observational lines of evidence provide at least par- 
tial support for such a picture. Measurements of ICM 
metal-abundances and their gradients (e.g. Loewenstein 
& Mushotzky 1996) are likely explained by feedback from 
star formation. Most of the stars in cluster elliptical galax- 
ies are of a nearly uniform old age (Bower, Lucey & Ellis 
1992; Kuntschner 2000, but see Trager etal 2000) and late- 
type galaxies within clusters are typically gas-poor and 
forming stars at a reduced rate compared to field galax- 
ies (Cayette et al. 1994, Chamaraux, Balkowski & Gerard 
1980). Even in MS1054, a very distant cluster at z = 0.83, 
much of the visible galaxy mass is contained in red, early- 
type galaxies (Van Dokkum et al. 1999). An active period 
of star formation in at least some regions of the z ~ 3 uni- 
verse is inferred from the ultraluminous sources detected 
by SCUBA (Barger, Cowie & Richards 2000), and the pop- 
ulation of Lyman-break galaxies, the distant equivalent of 
the local normal population, is forming stars at a rate sub- 
stantially higher than today (e.g., Steidel et al. 1999). 

In addition to winds driven by star formation, the jets 
of radio galaxies present another potential source of ICM 
heating. Chandra observations of low emission measure 
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bubbles in the ICM, aligned with the radio halos of cen- 
tral cluster sources (McNamara et al. 2000), implicate ra- 
dio jets as another source of additional ICM entropy. The 
fact that the radio galaxies appear more frequently in clus- 
ters at higher redshift (Owen 2000) suggests that heating 
from this population might be more important in the past. 

Of course, these and other observations also indicate 
that some late star formation and feedback does occur in 
cluster environments and this fact ultimately limits the ac- 
curacy and applicability of the approach used here. Also, 
the statistics of apparent gaps in the Ly-a forest of QSO 
spectra can be used to place limits on the volume fraction 
and epoch of heating to temperatures well above 10 4 K 
(Theuns, Mo & Schaye 2000). 

A principal benefit of a model where clusters form 
through the infall of preheated gas is that the physi- 
cal treatment of ICM evolution through gravitationally- 
driven shock heating is familiar territory for cosmological 
gas dynamic codes. Comparing the gravitational evolu- 
tion of a single X-ray cluster derived from twelve nearly- 
independent codes, Frenk et al. (1999) find agreement in 
structural properties at characteristic levels of ~ 5% for 
dark matter and ~10% for the ICM. 

The price for operating in a regime where the numerical 
accuracy has been calibrated is an approximate physical 
treatment. Neglect of radiative cooling is a particular con- 
cern, but the effect of radiative cooling and subsequent star 
formation on a cluster's thermal history is a complex prob- 
lem. Observations indicate that ICM cooling occurs fre- 
quently in clusters, but the phenomenon is restricted to a 
core region containing^ 1% of the collapsed (virial) cluster 
mass. However, recent spectroscopic studies which resolve 
cooling flows (David et al. 2000, Tamura et al. 2000) indi- 
cate that the phenomenon is confined to an even smaller 
volume than previously thought. 

At early times, cooling is an important ingredient of the 
galaxy formation process (Rees & Ostriker 1977; White & 
Rees 1978). So important, in fact, that sources of heating 
must be introduced to stabilize its effects on small scales 
(White & Frenk 1991; Cole 1991; Blanchard, Valls-Gabaud 
& Mamon 1992). Attempts to model the full problem of 
the galaxy formation process within clusters are limited 
by poor knowledge of the correct parameterizations for 
star formation and feedback at arbitrary epochs. Semi- 
analytic methods (Kauffman, White & Guiderdoni 1993; 
Baugh, Cole & Frenk 1996; Wu, Fabian & Nulsen 1998, 
1999; Somerville & Primack 1999) are best able to explore 
the large parameter space associated with these processes, 
but the lack of a direct solution for the spatial distribution 
of material limits the ability of this approach to predict de- 
tailed structural observables. Bryan (2000) has presented 
a model which shows that the elimination of low-entropy 
gas reproduces the luminosity-temperature scaling rela- 
tion. 

An initial entropy excess prevents gas from falling into 
the dark matter-dominated potential well to the extent 
it would in the purely gravitational infall model, thus re- 
ducing the central gas density. Core entropy 1 is related 
to temperature and density by S core ~ T/ p core 2 / 3 . If the 
temperature is set by the collisionless dark matter through 

1 In this paper, an equivalent entropy S is expressed in units keV cm 
where 7 = 5/3 for a monotonic gas (appropriate for a fully ionized pla: 
an ideal gas equation of state P = R p T. 



the virial theorem T oc M^o (M 200 is the mass enclosed 
in a sphere defined as containing a density contrast of 200 
with respect to the critical density), then the core density 
of initially preheated material will satisfy 

Pcore < M200 / Sf 2 . (1) 

Higher levels of preheating will produce lower core densi- 
ties and, for a given preheating amplitude Si, low mass 
clusters will have less dense cores than high mass clusters. 
This differential effect results in a steepening of the slopes 
of the relations between density-dependent cluster observ- 
ables (isophotal-size, luminosity and ICM mass) and tem- 
perature. 

Previous simulations have found preheating to pro- 
duce more realistic clusters in terms of the luminosity- 
temperature relation (e.g. Narvarro, Frenk & White 1995; 
Pierre, Bryan & Gastaud 1999). Metzler (1995) demon- 
strated a similar result using simulations that explicitly in- 
cluded the feedback of energy and mass from cluster galax- 
ies. We seek here to systematically investigate preheating's 
effect on the size-temperature, luminosity-temperature, 
and ICM mass-temperature relations. Although not en- 
tirely independent observables, the functional dependence 
on the ICM density varies for these measures. Each probes 
a differently-weighted moment of the radial density profile. 
There is no guarantee, therefore, that a single preheating 
level Si will provide a simultaneous match to all three re- 
lations. 

We use a total of 84 moderate resolution numerical ex- 
periments consisting of twelve initial configurations run at 
six different levels of preheating plus a set evolved without 
preheating for comparison. Since we lack a detailed mech- 
anism for the preheating, all simulations arc initiated with 
the gas held at a fixed, elevated temperature that corre- 
sponds to adiabatically evolved temperatures of ~ 10 6 K 
at z~3. 

In Section ||, the simulations are detailed. The effect 
of preheating on the major cluster relations are presented 
in Section Bl In Section |4| the range of entropies allowed 
by observation are determined and compared to other pre- 
heating studies. The effects on evolution, the virial mass- 
temperature relation, and gas fractions are explored in 
Section [|. 

2. TECHNIQUES 

2.1. The Simulations 

The simulations are run using the Lagrangian code 
P3MSPH (Evrard 1988). The parameters of the models 
we employ are similar to that used for the P3MSPH con- 
tribution to the Santa Barbara cluster comparison study 
(Frenk et al. 1999). Initial conditions for our clusters are 
produced using the same multi-step procedure used in the 
Frenk et al. study. First, purely gravitational N-body runs 
are used to identify an ensemble of clusters for resimulation 
with gas dynamics. The final resimulation uses a combi- 
nation of collisional and collisionless particles with a full 
gas dynamic treatment for a (Lagrangian) subset of the 
volume that comprises the cluster. An intermediate, low 

2 . The usual thermodynamic definition of entropy is s = cy ln(P/ p 7 ), 
3ma), relates to the entropy of this paper by S = e^ s ^ Cv ^ / R, assuming 
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resolution N-body model is used to identify the Lagrangian 
region to be treated with gas dynamics. 

The clusters are formed in a cold dark matter cosmol- 
ogy (Peebles 1982; Blumenthal et al. 1984) dominated 
by a non-zero cosmological constant, ACDM. The mod- 
els assume a flat (CI — 1) geometry with the following 
parameters: Q m = 0.3, Qa = 0.7, fib — 0.03, as — 1.0, 
T = 0.24 and h = 0.8. The Hubble constant is defined 
as 100 h km s^ 1 Mpc^ 1 ; and as is the power spectrum 
normalization on 8/i _1 Mpc scales. Clusters form from 
initial density perturbations which are Gaussian random 
fields consistent with a CDM transfer function specified by 
the shape parameter, T = fth (Bond & Efstathiou 1984). 
Initial conditions are constructed using Zcl'dovich's for- 
mulation (e.g. Efstathiou et al. 1985). 

The simulations begin with two 128 3 N-body runs which 
represent cubic regions 366 Mpc on a side from a CDM 
power spectrum. Clusters that form in this region have 
masses as large as 10 15 M© containing ~10 3 particles. We 
identify six clusters to resimulate from the final configura- 
tions of the two volumes: the two most massive clusters, 
a random pair less massive by a factor 3 from the mean of 
the first pair, and another random pair reduced in mass 
by a further factor 3. This results in a resimulated sample 
comprised of twelve clusters covering roughly a decade in 
mass. 

The density field in cubic regions centered on the chosen 
clusters' initial states are extracted from the original sim- 
ulation and placed onto a higher resolution 64 3 grid. The 
size of the regions ranges from 50-100 Mpc and scales as 
the cube root of the mass enclosed within the turn around 
radius of the resimulated cluster. High frequency modes 
of the density field not sampled by the original simulation 
are then added, up to the new limiting Nyquist frequency. 
From this set of initial conditions the cluster is evolved in 
a purely N-body simulation of 32 3 particles, using alter- 
nate sampling of the 64 3 density field. Particles in this 
simulation that lie within a density contrast of 5 with re- 
spect to the background are used to define a Lagrangian 
mask. For the final gas dynamic resimulation, masked 
locations in the 32 3 subsampled field are expanded by a 
factor of two in resolution, generating an effective 64 3 res- 
olution within the non-linear regions of the cluster. The 
high-resolution inner regions contain between 20,424 and 
26,064 dark matter and gas particles in equal numbers, 
where the gas particles are placed to trace the dark matter. 
The inner regions are surrounded by 29,510 to 30,215 low- 
resolution dark matter particles that are each eight times 
the mass of a combined high-resolution gas and dark mat- 
ter pair. This treatment allows the inclusion of both tidal 
effects from the surrounding large scale structure and the 
gas dynamics of the virial region in an economical way. 
The tidal particle contamination inside r2oo at z = is 
< 0.1% by mass. In addition, the cluster resimulations 
have mass resolutions which are similar fractions of the 
cluster virial masses despite spatial resolutions that range 
by a factor two, from 125-250 kpc. 

Preheating is achieved by specifying an initial hot tem- 
perature for the baryons. The initial epoch is defined by a 
linear growth factor of 16 to the present day, resulting in 
an initial redshift z t = 20.82 for ACDM. The initial tem- 
perature of the gas in the six preheating models ranges 
from 1.5 — 9.5xl0 7 K compared to 10 4 K in infall models. 



Since we expect that this heating actually occurs at z ~ 3, 
this method of preheating is valid only if a small amount of 
gas has collapsed at the true epoch of preheating. In these 
12 clusters 4 - 19% of the gas particles which eventually 
collapse have done so at z ~ 3. 

Table 1 lists the temperatures and corresponding initial 
entropies Si for the seven models. The latter are derived 
using the cosmic mean baryon density at the initial epoch. 
Since the simulations have density perturbations at the 
outset, the entropy floor is not perfectly flat. Also listed 
in Table 1 are values of the gas temperature at z = 3 that 
would place gas at the background density on the same 
initial adiabat. Energetically, at this epoch, the tempera- 
tures correspond to values between 0.07 and 0.4 keV per 
particle. 

Table || lists basic properties of the complete set of 84 
clusters at z — 0. Groups are numbered by mass beginning 
with the highest mass group from each of the two initial 
N-body models (labeled "a" and "b"). Individual runs are 
labeled by combining their model and group numbers with 
the initial entropy level (e.g. , al90S3). 



TABLE 1 
Model Properties 



Model 


Si 


T a diab{z = 3.0) 


Ti{z = 20.82) 




[keV cm 2 ] 


[K] 


[K] 


S6 


335.4 


5.3 x 10 B 


9.5 x 10 Y 


So 


221.8 


3.4 x 10 s 


6.0 x 10 7 


SI 


141.2 


2.2 x 10 6 


4.0 x 10 7 


S3 


105.9 


1.7 x 10 6 


3.0 x 10 7 


S2 


88.3 


1.4 x 10 s 


2.5 x 10 7 


SI 


53.0 


8.4 x 10 5 


1.5 x 10 7 


SO 


0.035 


5.6 x 10 2 


1.0 x 10 4 



2.2. Data Processing 

Mathiesen & Evrard (2001; hereafter ME01) find that 
there are significant differences (up to 20%) between the 
temperatures that theorists generally report (e.g. mass- 
weighted temperature) and their observed counterparts 
(which are based on plasma models). The values reported 
here build on the work of ME01 in order to be as close 
to observed values as possible. Luminosities and isopho- 
tal sizes are determined using the XSPEC mekal spec- 
tral emission model with 0.3 solar metallicity on a band- 
pass from 0.5 to 2.0 keV where the clusters are imaged 
at z=0.06. Unless otherwise stated, cluster temperatures 
used below are spectral values determined by the best fit 
of the cluster emission for a fixed photon count within 
rsoo to an isothermal mekal spectrum in a 0.5 to 10.0 keV 
bandpass. 

3. THE SCALING RELATIONS 

Clusters exhibit a number of scaling relations that 
lend insight into their physical nature. A cluster's 
X-ray isophotal size is tightly correlated with its emis- 
sion weighted mean temperature; the so-called size- 
temperature relation (i?/-T; Mohr & Evrard 1997; Mohr 
et al. 2000; hereafter MOO). The relationship between 
a cluster's luminosity and its temperature (Lx~T) has 
higher scatter, (David et al. 1993), reflecting its depen- 
dence on the ICM density distribution (the luminosity) 
and the total binding mass (temperature). This relation 
is important in determining the temperature function of 
a flux-limited sample of clusters. Because variations in 
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ICM mass fraction are rather modest (Mohr, Mathiesen & 
Evrard 1999; hereafter MME99), the ICM mass correlates 
with cluster temperature leading to an Mjcm~T relation. 
In addition to being useful in constraining cosmological pa- 
rameters, these relations can be used to constrain models 
of the interaction between galaxies and the ICM. 

3.1. Scaling Formalism 

The self-similar model of Kaiser (1986) is a useful start- 
ing point in developing the expected behavior of the scal- 
ing relations to be considered in this section. The model 
assumes a smooth, spherically symmetric distribution of 
gas and total mass about its center. We follow notation 
similar to that of Arnaud & Evrard (1999). Let M Aa be 
the total mass contained in a sphere (of radius r Aa about 
the center) that encompasses a mean density A c p c , where 
p c (z)=3H(z) 2 /8ttG is the critical density and H(z) the 
Hubble parameter of the universe at epoch z. We write 
the ICM density picm(t) in terms of the natural radial 
variable y=r/r& c 



PiCM(yrA c ) = fiCM A c p c g(y) 



(2) 



where explicit use of the ICM mass fraction within a den- 
sity contrast A c 



ICM 



M IC m(< ta c ) 
M Ac 



(3) 



sets the normalization of the structure function 
Ztfdyy 2 g(y) = l. 

A population of clusters that is strictly self-similar will 
have constant fjcM and a single, specific function g(y). 
Given the random and fully three-dimensional nature of 
cluster formation dynamics, strict self-similarity is an un- 
realistic expectation (Jing & Suto 1998; Thomas et al. 
2000). A more reasonable expectation is that clusters are 
a regular population whose gas fraction values and struc- 
ture functions will exhibit deviations depending, for ex- 
ample, on dynamical history. Searching for trends in the 
population with temperature, it is useful to employ mean 
gas fractions and structure functions 



fiCM(T) = (/jcm)o 
g(y\T) ee (g(y)) T 



(4) 
(5) 



where ( )t denotes an ensemble average over clusters with 
temperature T. 

Assume that the virial theorem holds, so that the mass 
Ma c and temperature are linked by 



M Ac = C m T° 



(6) 



where the exponent a m is 3/2 if spectral temperature is an 
unbiased measure of the mass-weighted value. Depending 
on the applied band pass, simulations display a bias that 
increases a m by ~ 10 — 20% (Mathiesen & Evrard 2001; 
and see §5.2 below). 

Finally, approximating cluster radial temperature pro- 
files as isothermal, we can write the temperature scaling 
relations for Mjcm, Lx and Ri that are expected from a 
regular cluster population described by mean gas fractions 



fiCM(T) and structure functions g(y\T). The ICM mass 
will scale as 

Micm = C m f ICM (T) T am (7) 

and the bolometric X-ray luminosity as 

Lx = C x ficM(T) 2 Q L (T)A(T)T a ™. (8) 

where A(T) is a dimensionless emissivity (see Appendix) 
and Ql(T) = 3 J q dyy 2 g 2 (y\T). The size at fixed isophote 
requires an additional assumption about the projected ra- 
dial profile of the emission. We assume a standard /3-model 
for which I(R)(xR 1 ~ 613 outside the core. For this case, the 
expected scaling of isophotal size is 

Ri = C R [f IC M 2 {T)Q I (T)UT)T 2a ^] 1 /W-V (9) 



with Qi{T)= j/ 1 ^ dr, h 2 (^e + rf) and 0.1 > £ 2 < 1. 
For the typical value (3 = 2/3 (Jones & Forman 1984; 
MME99), this becomes 



Ri 



Cr Jicm V \t) Q} /3 (T) A x / 3 (T) T 4q ™/ 9 . (10) 



For self-similar clusters and mass-weighted tempera- 
tures, familiar scalings Mjcm ocT 3//2 and LxotT 2 emerge 
for the ICM mass and luminosity. The isophotal size in this 
case scales as Ri cxT 2 / 3 . Equations (f?J), (0) and (||) demon- 
strate that deviations of observed slopes from these expec- 
tations can arise from a number of sources. Temperature- 
dependent mean gas fractions and structure functions may 
be largely responsible, but metallicity trends that affect 
A(T) and a virial slope a m ^ 3/2 will also tilt these rela- 
tions. 



3.2. The M 



ICM 



-T Relation 



Two-component /3-model analysis of 45 clusters in the 
X-ray flux-limited Edge sample leads to an observed 
Micm-T relation of the form (MME99) 

log(M ICM ) = (1.98±0.18)log(T 6 )+(13A2±0.03)-~log(h) 

(11) 

where Mjcm is the ICM mass within a density contrast 
of 500 in units of h~ l M Q and T 6 = T/6keV. The slope 
1.98 ± 0.18 is significantly different from the value of 3/2 
expected from a strictly self-similar population. 

Figure |l| displays the Mjcm~T relations at A c = 500 
for the simulated cluster sets at zero redshift. In this and 
subsequent figures of this section, the temperature is the 
spectral fit to emission within r^QQ. This radius is chosen 
because it is observationally accessible, but our results are 
fairly insensitive to this choice. Slopes differ at the level 
of 0.6 — 6% between r5oo and r2oo- The dotted line in each 
panel is the MME99 observational result scaled to h = 0.8. 
Increased levels of preheating lead to steeper slopes in the 
Micm~T relation; the six preheated models sweep right 
through the observed slope (note: S6 is not displayed to 
conserve space). Fits for the model sets are listed in Ta- 
ble ||. These fits are made only for clusters with kT s > 2.0 
keV, to reflect the range over which the MME99 fit was 
made. 



BIALEK, EVRARD & MOHR 



03 

o 

II 



10 



10 



14 



13 _ 



10 1 



10 



13 




2 3 4 5 6 789 2 3 4 5 6 789 2 3 4 5 6 789 

kT s [keV] 

Fig. 1. — The ICM mass— temperature relation at A c = 500 for 12 clusters run at 6 initial temperatures. The solid line 
is the best fit to clusters with kT s > 2 keV. Dotted lines in each panel show the MME99 observational result. 



Lower mass clusters have a larger core relative to their 
virial radius with respect to high mass clusters. This ra- 
tio of core radius to virial radius increases with increasing 
initial entropy. Consequently, the fraction of gas on adia- 
bats higher than the default self-similar case also increases. 
This has the effect of lowering the gas mass fraction within 
the virial region; implying a decreased Mjcm- 

This argument can also be cast in terms of the extropy 
floor discussed earlier. The entropy floor of the cluster ef- 
fectively 'plugs', in the sense of making Jeans stable, some 
amount of core material. Later infalling matter encounters 
this material at progressively larger radii as Si is increased, 
leading to spillover at larger radii. Lower mass clusters are 
affected by this to a larger degree because their virial adi- 
abat is lower than that of rich clusters, eq. (|l|). The fit for 
each set of models is listed in Table @. These fits are made 
only for clusters with kT s > 2.0 keV, to reflect the range 
over which the MME99 fit was made. 

Models S1-S4 have slopes in statistical agreement with 
the observed value. The zero-points are modestly higher 
than observed, which may be an indication that the global 
baryon fraction of 0.1 used in the simulations is~20% too 
high. Care must be taken in making this comparison be- 
cause the observed and simulated values scale differently 
with Hubble's constant (h~ 5 ^ 2 and ft -1 , respectively). Re- 
scaling to h = 0.7 would raise the observed values 22% 
higher than the simulated measures, essentially eliminat- 
ing the offset for model S3, for example. Since other factors 
like ICM dumpiness due to mergers (Mathiescn, Evrard 
& Mohr 1999) or a multiphase medium (Gunn & Thomas 
1996; Nagai, Sulkanen & Evrard 2000) can affect the ICM 
mass at similar (10 — 20%) levels, attempting to pin down 
the zero-point of this relation to a precision better than 
this remains problematic. 



3.3. The Lx — T Relation 

The Lx~T relation has received considerable attention 
from observers (e.g. Edge & Stewart 1991; David et al. 
1993) and computational model builders (Metzler 1995; 
Navarro, Frenk & White 1995; Allen & Fabian 1998; Eke, 
Navarro & Frenk 1998) over the last decade. An important 
contribution was made by Fabian et al. (1994; cf. Allen & 
Fabian 1998), who showed that excess core emission asso- 
ciated with cooling flows is the primary source of the large 
scatter observed in the Lx~T relation. The scatter is sig- 
nificantly reduced in analyses that either excise the core 
cooling flow regions (Markevitch 1998) or examine samples 
of clusters defined to possess weak cooling cores (Arnaud 
& Evrard 1999). 

Since our simulations do not include the effects of cool- 
ing, we will compare them to the results of the stud- 
ies that avoid cooling cores. Analyzing 26 clusters with 
accurate temperatures and inferred cooling rates M < 
100 M Q yr" 1 , Arnaud & Evrard (1999) find the corre- 
lation between bolometric luminosity and temperature 
L ocT q with q = 2.88 ± 0.15. Markevitch (1998) removes 
photons from the central 100 ft. -1 kpc when calculating 
both the temperature and luminosity and this produces 
q = 2.64 ± 0.27. We use an unweighted average of these 
two observational results, leading to an observed relation 
with intermediate slope 

log(Lx) = (2.76 ±0.15) log(T 6 ) + (44.48 ± 0.03) - 2log(h). 

(12) 

The bolometric Lx~T relations for the model sets at 
a redshift of zero are displayed in Figure 0. Power-law 
fits are listed in Table 0. Infall models (SO) exhibit a 
slope 1.44 ± 0.26 that is marginally inconsiste nt w ith the 
L (x T 2 self-similar expectation (see Section 3.1). The 
slope is biased by the two lowest temperature systems; 
each lies a factor of roughly three above a fit performed 
under their exclusion. Examination of the pair's density 
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Fig. 2. — The final bolometric X— ray luminosity vs. spectral temperature for the 12 clusters evolved from 6 initial 
temperatures. Temperatures are estimated within rsoo, luminosities are shown within r2oo (open squares) or within that 
radius but with the core emission (within 0.13r2oo) removed (crosses). Solid and dashed lines are fits to these respective 
data. The observational result for h = 0.8 is superimposed with a dotted line in each panel. 




2 3 4 5 6 789 2 3 4 5 6 789 2 3 4 5 6 789 

kT s [keV] 

Fig. 3. — The size-temperature relations for the ensemble are shown for an isophote level / = 3.0 X 10~ 14 j 
ergs/s/cm /arcmin in the 0.5 - 2.0 keV band. The solid line is the best fit to clusters with kT 3 > 2 keV. The ob- 
served relation is represented by the dotted line. 
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Table 3 

Fits" to the M IC m-T , L x -T & R r -T Relations at z = (h = 0.8) 





Micm-T 


Lx-T (core 


— included) 


Lx-T (core 


— extracted) 


Rr 


-T 


Model 


m 


b 


m 


b 


m 


b 


m 


b 


S6 


2.67(0.13) 


13.58(0.03) 


3.96(0.26) 


44.36(0.07) 


3.95(0.24) 


44.32(0.07) 






S5 


2.37(0.10) 


13.72(0.03) 


3.52(0.25) 


44.67(0.07) 


3.47(0.21) 


44.60(0.06) 


1.47(0.04) 


-0.36(0.01) 


S4 


2.07(0.08) 


13.80(0.02) 


3.03(0.24) 


44.84(0.07) 


3.00(0.18) 


44.75(0.05) 


1.20(0.07) 


-0.28(0.01) 


S3 


1.92(0.08) 


13.83(0.02) 


2.67(0.25) 


44.89(0.07) 


2.71(0.19) 


44.81(0.06) 


1.25(0.14) 


-0.24(0.02) 


S2 


1.86(0.07) 


13.85(0.02) 


2.55(0.24) 


44.93(0.07) 


2.57(0.17) 


44.84(0.05) 


1.15(0.10) 


-0.23(0.02) 


SI 


1.73(0.09) 


13.88(0.03) 


2.20(0.28) 


45.00(0.09) 


2.32(0.18) 


44.89(0.06) 


0.98(0.17) 


-0.18(0.02) 


SO 


1.69(0.12) 


13.95(0.04) 


1.56(0.38) 


45.09(0.14) 


2.02(0.13) 


44.96(0.05) 


0.74(0.05) 


-0.19(0.02) 


Observed 


1.98(0.18) 


13.66(0.03) 


2.76(0.15) 


44.67(0.03) 






1.02(0.11) 


-0.35(0.02) 






a Using the form log(X) 


= mlog(T s /6keV) + 6. 









profiles revealed that their core internal gas density profiles 
are steeper compared to the rest of the sample. A possible 
explanation of these enhancements is that they are tran- 
sient effects of fortuitously observed mergers (Roettiger et 
al. 1996). Exclusion of core emission, where the core is 
defined as a circular area of radius 0.13r2oo (Neumann & 
Arnaud 1999), yields a best fit L X ~T slope 2.04 ± 0.16, 
a value consistent with the analytic scaling. In Figure |[ 
crosses show the core-extracted Lx~T relation and the 
dashed line shows the best fit. 

Low-mass clusters, with lower virial entropy, feel the 
effects of preheating to a greater degree than high-mass 
clusters, resulting in a steepening Lx~T relation as Si is 
increased. The slopes of models S2-S4 are in agreement 
with observations. In the preheated models, no signifi- 
cant change in slope is seen by neglecting the core (see 
Table ||). The change in behavior compared to the SO 
models is expected when one considers the fact that the 
core density is bounded from above by the imposed initial 
entropy, eq. (@). 

The ability to match the Lx~T slope is not unique to 
the form of preheating assumed here. Metzler's (1995) 
models employing continuous energy feedback produce a 
slope of 2.96±0.05. The Lx~T relation is also explained 
by semi-analytic models such as the shock model of Cav- 
aliere, Menzi & Tozzi (1999) and models based on energy 
input and hydrostatic arguments (Wu, Fabian & Nulsen 
1998, 1999; Bower et al. 2000). Bryan (2000) matches 
the Lx~T relation with a model that removes low entropy 
ICM gas and assumes it to be trapped within galaxies. 

3.4. The Ri~T Relation 

The isophotal size Ri of a nearby cluster is tightly corre- 
lated to its temperature in a power law relation with slope 
near unity (Mohr et al. 2000, hereafter MOO) 

logRj = (1.02±0.11) log (T 6 ) -(0.44±0. 01)- log (h) (13) 

with Ri in units of Mpc. Mohr & Evrard (1997) showed 
that the observed slope is steeper than the range 0.61 < 
m < 0.81 found in sets of purely gravitational simula- 
tions of clusters in four different CDM cosmologies. The 
isophotal sizes from an ensemble of models experiencing 
continuous feedback showed a slope 0.99, consistent with 
observations. 

Figure || illustrates the effect of preheating on the Ri-T 
relation at zero redshift. The size Ri is determined at an 



isophote level / = 3.0 x 10~ 14 ergs/s/cm 2 /arcmin 2 in the 
0.5 - 2.0 keV band. In each panel, an individual cluster 
contributes three data points, one for each orthogonal pro- 
jection. The conversion from count rate to physical flux 
units is made by using PROS, assuming a Raymond-Smith 
spectrum with the published mean temperature Tx and -| 
solar abundances (Mushotzky & Loewenstein 1997). The 
angular diameter distance is used to convert angular size 
into physical size. I is also corrected for cosmological dim- 
ming, (1 + z) 4 . Table || lists the best fit power law param- 
eters for the simulations. Only clusters with kT s > 2 keV 
are used in the fit, so as to match the range of the observed 
dataset used by MOO. The fit is the average of 10,000 sets 
of 11 randomly drawn simulation images. Assuming that 
2 out of 3 projections can be treated as independent, each 
set excludes a random projection. 

Preheating leads to steeper slopes; the six models pro- 
duce a range of slopes that encompass the observed value. 
Model SI is in best agreement with observation, but mod- 
els S2 through S4 are also consistent. 

4. CONSTRAINING THE PREHEATING LEVEL 

Preheating succeeds in reproducing the slopes of the 
three observed scaling relations considered in the preced- 
ing section. Figure || summarizes the results by displaying 
the slope fit to each of the preheated sets along with the ob- 
servational constraints for the Micm~T , Lx~T and Rj—T 
relations. Simultaneous agreement in all three slopes can 
be satisfied within a range of initial entropy values near 
-lOOkeVcm 2 . 

We use a chi-squared measure based on the slopes of the 
three relations Ax 2 = J2i=i { m sira,i - m obSti ) 2 / (a 2 imi + 
a obs i) ~ Xmin f° r eacn set 01 models in order to constrain 
the range of allowed preheating. Figure 5] shows the chi- 
squared for each relation individually and their sum. Hor- 
izontal lines indicate the 68.3%, 90% and 99% confidence 
limits for one degree of freedom. From this analysis, the 
predicted level of preheating required to reproduce obser- 
vation is between 55-150 keV cm 2 at 90% confidence. 

This range overlaps values derived independently by 
Lloyd-Davies, Ponman & Cannon (2000). In an analysis 
of X-ray emission profiles of groups and clusters, they in- 
ferred a trend of decreasing central density with decreasing 
temperature, consistent with the constraint of a minimum 
central entropy in the range 34-82 keV cm 2 with 90% 
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Fig. 4. — The Micm~T , L x -T , and slopes as a func- 

tion of initial entropy. The allowed slopes based on observation are 
indicated by a pair of dotted lines. All errors are la. The open 
symbols in the Lx-T panel represent the core-extracted fits (the 
two symbols have been offset about the Sj value for clarity). 

certainty. Pierre, Bryan, and Gastaud (f999) produced a 
suite of simulated clusters at a preheated entropy of 141 
keV cm 2 and achieve results that are in agreement with 
the observed Lx~T relation (see also Pen 1999). Balogh, 
Babul & Patton claim success with a model in which the 
gas at the turnaround density is heated to 1.8 x 10 6 K at 
z = l, translating (for £l m = 0.3, h — 0.8) into an equivalent 
entropy value 170 keV cm 2 . 

As previously mentioned, preheating is not unique in its 
ability to solve the Lx~T scaling relation problem. Semi- 
analytical models (Wu, Fabian and Nulson 1999a, b; Bower 
et al. 2000; Tozzi & Norman 2000) typically add energy 
into the ICM in a more gradual fashion. Since the evolu- 
tionary history of the ICM would differ in different scenar- 
ios, the behavior of the distant cluster population should 
decide between ICM histories based on rapid or more grad- 
ual heating. We turn next to predictions for ICM evolution 
at z = 0.5 and z = l. 

5. EVOLUTIONARY PREDICTIONS 

Preheating affects the evolution of the ICM in a char- 
acteristic way. Observations of the distant cluster popu- 
lation can therefore be used to falsify the model or de- 
fine the limits of its applicability. Figures || through || 
give a representation of the effect of preheating on cluster 
al. The projected X-ray emissivity, thermal SZ signal, 
and emission-weighted temperature at redshifts 0, 0.5, 1 
and 2 are shown for the same cluster under three treat- 
ments: no heating (model SO), a successful level of pre- 
heating (S3) and an excessively heated case (S6). The 
surface-brightness and temperature maps are made using 
the mekal emission model and the images are sized to be 
2?*2oo on a side. The excessive model S6 shows the power 
of preheating to deplete the central ICM density of a clus- 
ter and severely limit its X-ray visibility at high redshift. 
The impact on the preferred model S3 is less dramatic. In 




40 50 60 70 809CL00 200 300 

Sj [keV cm 2 ] 

Fig. 5. — The A\ 2 measure for the scaling relation slopes is 
shown as a function of initial entropy. The dotted, heavy dashed, 
dot-dashed lines correspond to the Micm~T , Lx—T and Rj—T 
relations, respectively. The solid line indicates the \ 2 when consid- 
ering all three relations. The light dashed and solid lines correspond 
to the use of the core-extracted Lx—T fits. 

this section, we present the redshift evolution of the locally 
successful set of S3 models only. This model is chosen be- 
cause it lies mid-range in the allowed initial entropies and 
because it is the best fit for two of the three relations. 

5.1. Evolution of Scaling Relations 

Figures ^|, |l0| and [ll] show evolution in the Micm~T, 
Lx~T and Rj-T relations, respectively, for model S3 at 
redshifts 0, 0.5 and 1. Table [| provides fits to power 
law relations. The ICM mass and luminosity are mea- 
sured within r2oo at each redshift and the physical isopho- 
tal size is determined at a fixed intrinsic isophote level of 
I = 3.0 x 10" 14 ergs/s/cm 2 /arcmin 2 in the 0.5 - 2.0 keV 
band. 

The simulations predict little evolution in the Rj-T and 
Lx~T relations to z = 0.5. However, at z = 1.0, both rela- 
tions become ~ 30% steeper. Current observations of high 
redshift clusters appear consistent with these predictions. 
Mohr et al. (2000) find no significant evolution in their in- 
termediate redshift study (z w 0.2-0.5) of CNOC clusters 
imaged by ROSAT. There is no significant evolution in 
the Lx~T relation to redshifts ^0.5 (Mushotzky & Scharf 
1997; Fairley et al. 2000). Higher quality data from the 
coming Chandra and XMM/Newton archives will tighten 
the constraints on the degree of allowed evolution. 

Figure |9| illustrates that the Mjc m~T relation evolves 
toward lower ICM masses at higher redshift for fixed tem- 
perature clusters, as expected in the canonical evolution 
model. As with the Rj-T and Lx~T relations, no signifi- 
cant evolution in slope is expected at redshift 0.5 while a 
moderate steepening to Mjcm oc T' 2 ' 4±0 - 2 is seen at z=l. 
Because ICM masses within r2oo require very deep imag- 
ing, less observational information is available for this re- 
lation at high redshifts. The expectations of model S3 ap- 
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Fig. 9. — The ICM mass-temperature relation for the clusters of 
preheated model S3 at three redshifts. 

ppear consistent with observations (e.g. , Matsumoto et 
al. 2000). 

5.2. Virial Mass & ICM Mass Fractions 

The fact that the Mjcm~T relation of set S3 is steeper 
than the canonical 3/2 expectation arises from two sepa- 
rate factors. One is the relation between total mass and 
temperature and the other is the ICM mass fraction within 

T200- 

Non-preheated models follow a relation between total 
mass Ma c and the mass-weighted temperature T m that 
is consistent with expectations from the virial theorem 
H(z)M Ac oc Til 2 (e.g. Bryan & Norman 1998; MOO; 
ME01). When a 2 — 10 keV spectral temperature is em- 
ployed in place of the mass- weighted value, ME01 find that 
the relation steepens slightly to H(z)Ma c oc T} n Q . Our 
non-preheated set of models (SO) is consistent with this 
steepening. 

For the preheated set S3, Figure |l^ presents the rela- 
tion between total mass within r2oo versus spectral tem- 
perature at z = 0, 0.5 and 1. Power-law fits to the data, 
listed in Table 5, indicate that slopes of 1.66 ± 0.12 and 
1.73 ± 0.11 at z = and 0.5, respectively, are consistent 
with the spectral relation of the non- preheated models. At 
z = l, the slope steepens to 1.91 ± 0.11. Note that the in- 
tercept at T s — 6 keV remains stable, varying by only a 
few percent over the redshift range probed. 

Observational determination of this relation can be at- 
tempted using weak lensing masses. Current data are 
too noisy to discriminate between a slope of 1.5 and 1.6 
(Hjorth, Oukbir & Van Kampen 1999). Weak lensing 
masses are likely biased by the structures in which the 
clusters are embedded (Metzler et al. 1999); In addition 
to lensing masses, estimates based on galaxy kinematics 
(Girardi et al. 1998; Horner, Mushotzky & Scharf 1999) or 
the hydrostatic assumption with measured teni- 
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Fig. 10. — The luminosity— temperature evolution for the clusters 
of model S3. 




Fig. 11. 
0.5 and 1. 



3 4 5 

kT s [keV] 

The size-temperature relation for model S3 at z = 0, 



perature profiles are possible. Nevalainen, Markevitch & 
Forman (2000) follow the latter approach for a sample of 
nine galaxies, groups and clusters and find a relation with 
slope 1.79 ± 0.14, steeper by 3.6<r than the canonical 1.5 
expectation but more in line with the relation from set 
S3. Finoguenov, Rciprich & Boehringer (2000) analyze a 
larger sample and find a slope of l.78t^\og, consistent with 
the Nevalainen et al (2000) results and inconsistent with 
self similar expectation at more than 3a. 

The ICM mass-temperature relations are steeper than 
that for the total mass at all redshifts for set S3, because 
the ICM mass fraction decreases at lower cluster temper- 
atures. Figure Q_3| shows that the baryon fraction within 
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^200 decreases from 9% at high temperatures to about 6% 
at 2 keV. Note that all values are consistently lower than 
the cosmic baryon fraction of 10%, an effect generally seen 
even in non-preheated simulations (MME99; Frenk et al. 
1999). The fall-off at low temperatures is evidence of pre- 
heating's stronger impact on shallower gravitational po- 
tentials. 

Because cluster potentials tend to deepen with time, clus- 
ters evolve along this relationship toward values of higher 
ICM mass fraction. 

Preheating and its effects on the ICM mass fraction and 
the Mjcm-T relation is of direct relevance to the observ- 
ability of high rcdshift clusters based on their X-ray emis- 
sivity or the thermal Sunyaev-Zel'dovich effect. The latter 
is particularly well suited for high redshift cluster detec- 
tion. Using non-preheated simulations which with con- 
stant ICM mass fraction, Holder et al. (2000) show that 
proposed interferometric arrays would be capable of de- 
tecting clusters above a total mass limit of 10 14 ft. _1 M 
at essentially any redshift. Haiman, Mohr & Holder (2000) 
emphasize the power of such surveys to constrain cosmo- 
logical parameters f^M, and the equation of the state 
parameter w of the dark energy component. An analysis of 
the effects of preheating on these yields and possible cos- 
mological constraints is underway (Mohr et al. in prep.). 

6. CONCLUSIONS 

We have used a set of 84 moderate resolution gas dy- 
namic simulations of cluster evolution to systematically 
investigate the effects of preheating on the local Ri~T , 
Lx~T and Mjcm~T relations. We use a plasma emis- 
sion model to estimate luminosities and determine spec- 
tral temperatures from fits to the overall cluster emission. 
We find that preheated models with initial entropy in the 
range 55 — 150 keV cm 2 reproduce the slopes of the ob- 
served scaling relations. 

There is agreement between the observations of Lloyd- 
Davies et al. (2000) and these simulations in the allowed 
range of initial entropy. Although the range is currently 
generous, the new generation of X-ray satellites should 
provide information to further constrain the allowed range 
or potentially falsify the model. Better understanding of 
the ICM entropy aids in the deduction of its history and 
places restrictions on the mechanisms that may be respon- 
sible for preheating. 

Examining evolutionary effects in a set with Si = 
106 keV cm 2 , the models predict modest steepening in the 
three scaling relations at z = 1.0. The predictions are 
currently in agreement with 



Table 5 

Fits™ to h(z)M2oo vs. spectral temperature for set S3 



z 


m 


b 


0.0 


1.66(0.12) 


14.98(0.04) 


0.5 


1.73(0.11) 


14.95(0.04) 


1.0 


1.91(0.11) 


14.96(0.04) 



a \og(h{z)M 200 ) = mlog(T 3 /6kcV) + fe. 
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Fig. 12. — Redshift-corrected total mass within r200 vs. spectral 
temperature for clusters in set S3. 
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Fig. 13. — ICM mass as a fraction of total mass vs. spectral tem- 
perature inside r200 for clusters of model S3 at z = 0.0, 0.5 and 1.0. 
The global value is Qb/Qm =0.1. 

the observed lack of evolution in the Lx~T relation to 
z ~ 0.5. Deeper observations {e.g. , Stanford et al. 2000) 
from Chandra and XMM-Newton should help reveal the 
epoch where the scalings depart from the simulation's ex- 
pectations. 
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Table 4 

Fits" to the Ri-T , L x -T , Mjcm-T relations for preheating set S3 (h = 0.8) 





Micm-T 


Lx 


-T R r 


-T 


z 


m 


b 


m 


b m 


b 


0.0 
0.5 
1.0 


1.90(0.07) 
1.97(0.09) 
2.43(0.23) 


13.83(0.02) 
13.62(0.03) 
13.49(0.08) 


2.67(0.25) 
2.59(0.25) 
3.25(0.43) 


44.89(0.07) 1.16(0.12) 
44.81(0.08) 1.27(0.19) 
44.85(0.15) 1.49(0.19) 


-0.25(0.02) 
-0.27(0.05) 
-0.25(0.03) 






a Using the fo: 


rm log(X) = 


mlog(T s /6 keV) +b. 
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APPENDIX 

Let Ma c be the total mass contained in a sphere (of radius r/± c about the center) that encompasses a mean density 
A c p c , where p c (z)=3H (z) 2 / 8ttG is the critical density and H{z) the Hubble parameter of the universe at epoch z. Write 
the ICM density Picm{t") in terms of the natural radial variable y=r/r^ c , structure function g(y), and the ICM gas 
fraction within rA c , as in eqn's (||) and (0). 

Clusters have no well-defined edge, but their hydrostatic regions are reasonably well bounded by r2oo- We normalize 
the total cluster luminosity to be that interior to this radius 

L x = 4tt / drr 2 n e (r) n H (r) A(T(r)) (1) 
Jo 

where A(T(r)) is an appropriately normalized emissivity dependent only on temperature and n e = picm / Pe^np and 
tin = Picm I M-ff TO p are the number densities of free electrons and protons. Assuming that the ICM is isothermal, we can 
define a dimcnsionless emissivity A(T) = A(T)/A(10 keV) by arbitrarily normalizing to emission at T = 10] keV. Using 
equations (^|) and (|^), the luminosity can be rewritten as 

L x = C x fiCM 1 Ql A(T) T«™ (2) 

where Cx carries dimension and depends only on fundamental constants and Ql = 3 dyy 2 g 2 (y). Note that is 
equivalent to the clumping factor < p 2 ICM > / < Picm > 2 , where the angle brackets denote the volume average over 
the cluster atmosphere. It is thus a structure factor which depends solely on the gas density shape and characterizes the 
concentration of the gas distribution. 

Allowing expectations for ficM and Ql to vary with cluster temperature produces eq. (n) ■ To see the origin of the 
luminosity's canonical infall model temperature dependence, a set of additional assumptions are traditionally made: (i) 
pure bremsstrahlung emission (A(T) cx T 1 / 2 ), (ii) virial equilibrium (M200 oc T 3 ^ 2 ), (ni) structurely identical clusters 
(Ql(T) = Ci) and (iv) constant gas fraction (ficM — C2). This leads to the expected scaling relation between 
luminosity and temperature: Lx oc T 2 . 

The scaling for the Rj-T relation is outlined in Mohr et al. (2000). Assuming no emission beyond V2oo, then the surface 
brightness at r2oo is formally zero. Consider the surface brightness at radius £r2oo where £ is near, but below, unity 



1 rr2aa\/l-Z 2 

S(£r 2 oQ) = dl n e (r) n H (r)A(T) (3) 

2tt(i + zy j- rwo ^e 



where r = U (£r2oo) 2 + 1 2 - Rewriting using dimensionless variables as in the luminosity case above, we find 

£(£r 20 o) = Ci ficM 2 Qi A(T) T a ^ (4) 



where Cj is a constant and Qj = ^ dr/ g 2 (^/ rf + £ 2 ) . 

Assume that the j3- model (Cavaliere & Fusco-Femiano 1976) describes the surface brightness profile 




-3/3+i 

(5) 
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Outside the core, the surface brightness scales as a power law in radius E(i?) oc R- 6I3+1 . In this regime, the surface 
brightness at an isophotal distance, X(i?j), can be related to the surface brightness at radius ^200 



Rl = £ r 200 



1/(6/3-1) 



(6) 



This adds a constraint on the shape parameter Qjoc^ 1 6,3 . Using eq. (^J) and introducing temperature dependence in 
fiCM and Qi then leads to the result eq. (^|) of §3.1. For a typical observed value of (3 = 2/3 along with canonical 
assumptions for the other parameters, the predicted size temperature relation is Ri oc T 2 / 3 . 
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Table 2 

Cluster properties (A c = 200, z = 0) 



Model 


L a 


M tot b 


Mjcm c 

JO 1V1 


T 


^ d 
(7 dm 


r£ C 


T s l 


Lx bol^ 




alS6 


72.53 


8.48 


7.02 


0.828 


-85 


5.72 


6.24 


5.65 


0.3719 


a2S6 


72.53 


3.01 


1.88 


0.625 


-303 


3.06 


2.99 


0.408 


- 


a31S6 


53.34 


3.90 


2.54 


0.651 


-18 


3.69 


3.94 


0.937 


- 


a34S6 


42.86 


1.76 


0.872 


0.495 


-91 


2.23 


2.48 


0.164 


- 


al90S6 


37.89 


1.37 


0.395 


0.288 


-117 


1.89 


2.07 


0.0377 


- 


a260S6 


31.04 


0.896 


0.318 


0.355 


-30 


1.52 


1.80 


0.0498 


- 


blS6 


82.43 


10.4 


8.80 


0.846 


-311 


6.41 


5.96 


4.18 


- 


b2S6 


86.26 


11.7 


10.4 


0.889 


84 


7.81 


8.23 


8.13 


- 


b31S6 


51.56 


5.31 


3.31 


0.623 


-208 


4.43 


4.96 


2.64 


- 


b34S6 


48.23 


3.52 


1.82 


0.517 


-157 


3.35 


3.59 


0.476 


- 


b233S6 


40.91 


1.48 


0.485 


0.328 


-53 


2.06 


2.27 


0.0558 


- 


b255S6 


32.99 


1.30 


0.493 


0.379 


-67 


1.85 


2.16 


0.0692 


- 


alS5 


72.53 


8.96 


7.89 


0.881 


-51 


5.91 


6.24 


8.39 


0.4308 


a2S5 


72.53 


3.58 


2.66 


0.743 


-331 


3.10 


3.13 


0.889 


- 


a31S5 


53.34 


4.17 


3.35 


0.803 


-37 


3.65 


3.94 


1.93 


0.2215 


a34S5 


42.86 


1.95 


1.41 


0.723 


-80 


2.27 


2.48 


0.470 


- 


al90S5 


37.89 


1.43 


0.655 


0.458 


-61 


1.85 


2.16 


0.122 


- 


a260S5 


31.04 


0.952 


0.497 


0.522 


-91 


1.48 


1.72 


0.131 


- 


blS5 


82.43 


11.4 


10.2 


0.895 


-316 


6.23 


5.44 


5.69 


0.2519 


b2S5 


86.26 


12.1 


11.4 


0.942 


69 


7.65 


7.50 


10.2 


0.4759 


b31S5 


51.56 


5.38 


4.34 


0.807 


-22 


4.69 


5.19 


7.71 


0.3442 


b34S5 


48.23 


3.74 


2.73 


0.730 


-62 


3.58 


3.94 


1.58 


0.1922 


b233S5 


40.91 


1.58 


0.775 


0.491 


-87 


1.99 


2.27 


0.160 


- 


b255S5 


32.99 


1.41 


0.753 


0.534 


-52 


1.85 


2.16 


0.239 


- 


alS4 


72.53 


8.88 


8.08 


0.910 


-74 


5.79 


5.96 


9.70 


0.4395 


a2S4 


72.53 


3.70 


2.96 


0.800 


-343 


2.93 


2.99 


1.12 


- 


a31S4 


53.34 


4.34 


3.75 


0.864 


-63 


3.57 


3.94 


2.73 


0.2908 


a34S4 


42.86 


2.06 


1.76 


0.854 


-50 


2.36 


2.48 


0.822 


0.1435 


al90S4 


37.89 


1.50 


0.944 


0.629 


-64 


1.79 


2.07 


0.283 


- 


a260S4 


31.04 


1.23 


0.760 


0.618 


-188 


1.46 


1.72 


0.305 


0.0922 


blS4 


82.43 


11.5 


10.2 


0.887 


-309 


5.96 


5.19 


6.12 


0.3221 


b2S4 


86.26 


13.0 


12.1 


0.931 


130 


7.76 


7.50 


15.4 


0.5218 


b31S4 


51.56 


5.46 


4.93 


0.903 


30 


4.75 


4.96 


12.1 


0.3749 


b34S4 


48.23 


3.98 


3.29 


0.827 


-17 


3.70 


3.94 


2.61 


0.2854 


b233S4 


40.91 


1.65 


1.07 


0.648 


-86 


1.93 


2.16 


0.345 


- 


b255S4 


32.99 


1.62 


1.05 


0.648 


-87 


1.84 


2.07 


0.507 


0.1283 


alS3 


72.53 


8.88 


8.24 


0.928 


-53 


5.80 


5.96 


10.5 


0.4500 


a2S3 


72.53 


2.51 


3.50 


0.839 


-216 


3.35 


3.59 


2.06 


0.2398 


a31S3 


53.34 


2.72 


3.80 


0.892 


-40 


3.41 


3.59 


2.20 


0.3031 


a34S3 


42.86 


1.37 


1.91 


0.884 


-71 


3.40 


2.48 


1.08 


0.1873 


al90S3 


37.89 


0.853 


1.19 


0.735 


-71 


1.83 


2.07 


0.501 


0.0990 


a260S3 


31.04 


0.716 


1.00 


0.694 


-174 


1.59 


1.80 


0.598 


0.1485 


blS3 


82.43 


7.16 


10.0 


0.893 


-300 


5.81 


4.96 


6.37 


0.3608 


b2S3 


86.26 


8.70 


12.2 


0.938 


110 


7.81 


7.50 


16.1 


0.5418 


b31S3 


51.56 


3.74 


5.22 


0.944 


116 


4.69 


4.73 


12.8 


0.3905 


b34S3 


48.23 


2.48 


3.46 


0.865 


41 


3.76 


4.12 


3.49 


0.3139 


b233S3 


40.91 


0.922 


1.29 


0.741 


-42 


1.98 


2.27 


0.585 


0.1122 


b255S3 


32.99 


0.910 


1.27 


0.747 


-13 


1.94 


2.16 


0.869 


0.1825 


alS2 


72.53 


8.88 


8.32 


0.937 


-78 


5.71 


5.96 


10.7 


0.4495 


a2S2 


72.53 


3.75 


3.18 


0.848 


-369 


2.79 


2.72 


1.34 


0.0808 


a31S2 


53.34 


4.29 


3.80 


0.886 


-39 


3.43 


3.59 


3.15 


0.3145 


a34S2 


42.86 


2.05 


1.82 


0.888 


-68 


2.28 


2.37 


0.991 


0.1788 


al90S2 


37.89 


1.59 


1.26 


0.792 


-87 


1.82 


2.07 


0.593 


0.1315 


a260S2 


31.04 


1.50 


1.12 


0.747 


-158 


1.65 


1.88 


0.920 


0.1727 


blS2 


82.43 


11.1 


9.84 


0.886 


-293 


5.72 


4.96 


6.26 


0.3780 


b2S2 


86.26 


13.0 


12.2 


0.938 


104 


7.79 


7.50 


17.7 


0.5484 



Effects of Preheating in X-ray Clusters 



Table 2 — Continued 



Model 


L a 


M tot h 


Mjcm c 


T 


0~DM 


rp C 

■*■ m 




Lx,boi s 


Ri h 


UOIDZ 




o.uz 


o.ou 




1 OS 
luo 






19 4 

_L Z 




V.O /ICO 


/I Q OO 
40. zo 


a no 


0.00 


U.oOo 


00 
zo 


O. i i 


A 1 O 
4.1Z 


q 01 

o.ol 


n 0000 

U.ozoo 


KOQQCO 


a n oi 


1 70 
1. / z 


1 CO 
l.DO 


U.oyU 


00 
-00 


1 on 
i.yu 


1 ft 

z.10 


U.000 




DZOODZ 


on on 


1 7/1 

1. f 4 


1 07 


U. f f 





1 OO 

i.yy 


Olft 

z.io 


1.10 


n onnn 

u.zuuy 


albl 


72.53 


00 
o.oZ 


o o ,) 

8.24 


n no /i 
0.934 


^ *7 

-47 


C TO 

O.YO 


c ftn 
5.69 


101 
12.1 


n /I 0£? K 


a2bl 


TO C O 

/2.53 


<~t A A 

ZA4 


2.04 


n 00c 
0.836 


1 n£ 
-196 


01 
2.21 


07 
2.2 ( 


n 1 1 
0.811 


n n /i to 


a31bl 


53.34 


A A O 

4.42 


O OT 

3.87 


n o 7c 
0.876 


1 

-18 


cn 


c n 
3.59 


ftn 
3.60 


n a c 


a34bl 


,11 o£; 

42.00 


o 1 n 
2.11) 


1 n a 
1.94 


n no a 
0.924 


A 1 

-41 


O O A 

Z.Z4 


O O 7 

2.37 


1 C7 

1.57 


n oo£?n 


a!90bl 


07 on 
37.89 


1 71 

1.71 


1 Af\ 

1.49 


n 071 
0.8 ( 1 


c 
-53 


i.yu 


1 ft 
2.16 


1 no 
1.03 


n 01 on 
U.212U 


azoUbl 


31.04 


1.59 


1.36 


n o c c 
0.855 


-152 


1 TO 

1.72 


1 00 
1.88 


1 on 
1.80 


n onoo 


blSl 


OO A O 

82.43 


11.2 


9.60 


n cr 7 

0.857 


-257 


5.63 


/) nft 

4.96 


ft 7ft 

6.76 


n AO TO 

0.4378 


UOOI 


86.26 


13.4 


12.5 


a no 

0.933 


no 

92 


7 on 

7.80 


1-7 r n 

7.50 


on 

20.3 


n c n /i t 

0.5947 


b31bl 


51.56 


5.84 


5.48 


0.938 


143 


4.57 


4.32 


14.0 


0.4281 


b34Sl 


48.23 


3.98 


3.52 


0.884 


-81 


3.64 


4.12 


4.84 


0.3605 


b233Sl 


40.91 


1.74 


1.52 


0.874 


-7 


1.96 


2.16 


1.05 


0.2162 


b255Sl 


32.99 


1.78 


1.54 


0.865 


2 


2.01 


2.16 


1.81 


0.2389 


alSO 


72.53 


8.96 


8.32 


0.929 


-72 


5.55 


5.19 


15.4 


0.4547 


a2S0 


72.53 


3.41 


2.81 


0.824 


-308 


2.52 


2.48 


1.58 


0.2304 


a31S0 


53.34 


4.40 


3.84 


0.873 


-26 


3.28 


3.28 


4.56 


0.3349 


a34S0 


42.86 


2.11 


1.93 


0.915 


-73 


2.13 


2.07 


2.59 


0.2410 


al90S0 


37.89 


1.86 


1.66 


0.892 


5 


1.95 


2.07 


1.79 


0.2283 


a260S0 


31.04 


1.75 


1.55 


0.886 


-136 


1.53 


1.64 


5.05 


0.2358 


blSO 


82.43 


10.7 


9.12 


0.852 


-202 


5.44 


4.32 


7.33 


0.4339 


b2S0 


86.26 


13.6 


12.6 


0.926 


-77 


7.73 


7.17 


28.0 


0.5353 


b31S0 


51.56 


5.95 


5.41 


0.909 


147 


4.19 


3.42 


15.4 


0.3923 


b34S0 


48.23 


4.06 


3.58 


0.882 


-42 


3.59 


3.59 


6.78 


0.3337 


b233S0 


40.91 


1.80 


1.59 


0.883 


18 


1.94 


2.07 


1.94 


0.2373 


b255S0 


32.99 


1.92 


1.64 


0.854 


3 


1.85 


1.80 


3.84 


0.2228 



a [ h- 1 Mpc ] 
b [ h- 1 10 14 M Q ] 
c [ h- 1 1O 13 M ] 
d [ km/s ] 
C [K] 
f [K] 

s[ h 10 44 erg/s] 

h [ /i~ 4 Mpc ], dashed entries result from no emission above the chosen isophote 



BIALEK, EVRARD & MOHR 



15 



Fig. 6. — X-ray surface-brightness (SB) images for cluster al are shown for non-preheated (SO, left), moderately preheated (S3, middle) 
and excessively preheated (S6, right) cases at four redshifts (0.0, bottom; 0.5; 1.0; 2.0, top). 

Fig. 7. — X-ray temperature maps for cluster al are shown for non-preheated (SO, left), moderately preheated (S3, middle) and excessively 
preheated (S6, right) cases at four redshifts (0.0, bottom; 0.5; 1.0; 2.0, top). 

Fig. 8. — Sunyaev-Zcl'dovich effect (SZ) Compton Y images for cluster al arc shown for non-preheated (SO, left), moderately preheated 
(S3, middle) and excessively preheated (S6, right) cases at four redshifts (0.0, bottom; 0.5; 1.0; 2.0, top). 
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